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Abstract 

High-quality superconducting oscillators have been successfully used for quantum control and 
readout devices in conjunction with superconducting qubits. Also, it is well known that squeezed 
states can improve the accuracy of measurements to subquantum, or at least subthermal, levels. 
Here we show theoretically how to produce squeezed states of microwave radiation in a super- 
conducting oscillator with tunable parameters. The circuit impedance, and thus the resonance 
frequency, can be changed by controlling the state of an RF SQUID inductively coupled to the 
oscillator. By repeatedly shifting the resonance frequency between any two values, it is possible 
to produce squeezed and subthermal states of the electromagnetic field in the (0.1-10) GHz range, 
even when the relative frequency change is small. We propose experimental protocols for the veri- 
fication of squeezed state generation, and for their use to improve the readout fidelity when such 
oscillators serve as quantum transducers. 
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The problem of quantum measurements has recently attracted renewed attention. In 
quantum mechanics, the extraction of information from a quantum system produces an 
unavoidable disturbance on it. If the object is initially in an eigenstate of the measured 
observable, a quantum non-demolition (QND) measurement can be realized, where this 
disturbance is minimal [1]. 

One known type of detector for QND measurements is the so-called parametric trans- 

n 

ducer p|. A key element of a parametric transducer is an optical or radio-frequency auto- 
oscillator. A transducer coupled to a quantum system of interest is designed such that the 
behaviour of the quantum system changes the phase and/or amplitude of its oscillations. 
The phase (amplitude) shift provides information about the quantum system's dynamics. 
With the recent development of superconducting qubits , this approach was success- 
fully applied to their study. In particular, transducers (high-quality superconducting tank 
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circuits) were used to measure the state of superconducting flux qubits 

It is known that the noise of detectors can be decreased, even below the standard quantum 
noise level, by employing squeezed states 0, Q, Q] . In this paper we show that a supercon- 
ducting parametric transducer allows a natural application of this approach, since it can 
be used both to produce squeezed states and to use them in order to minimize quantum 
fluctuations. An immediate result of this method's application would be a way to suppress 
the effective noise temperature of the next-stage amplifier, at least to the nominal temper- 
ature of the cooling chamber. We emphasize that existing experimental techniques should 
be sufficient for the realization of our proposal. 

Squeezed states are quantum states in which the dispersion of one variable can 

be below the standard quantum limit (SQL). More specifically, if a system is described by 
a pair of dimensionless conjugate variables, Q and P, it is in a squeezed state if for some 
times 



{AQ') ^ {Q') - {Q)' < 1/2. (1) 

The uncertainty principle requires that (AP^)(A(5^) > 1/4, so when one variable is 
squeezed, the dispersion of the conjugate variable increases. Squeezed states were introduced 



in the context of quantum optics (see, e.g.. Re: 



number of other systems, including polaritons [9|,ll0||, phonons [ll|,ll4|, Josephson junctions 



13H and molecular oscillations 




), but were since investigated in a 



14| . Squeezed states have been successfully generated in 



a Josephson parametric amplifier 15i. A special interest in these states is due to their 

riQ 

usefulness in obtaining SQL resolution in imaging and measurement (see, e.g., [a, 17[). Their 
classical analog can be used to obtain subthermal resolution in mechanical measurements 

Q. 

As described in, e.g., Ref. 17|, if a harmonic oscillator is in a coherent state (i.e., a state 
with equal and minimal uncertainties: (AP^) = (AQ"^) = 1/2), a sudden change of the 
oscillator frequency would create a squeezed state (while an adiabatic change would not). 
The degree of squeezing is given by the ratio 

X = u/uq (2) 

of the oscillator frequencies before and after the shift, or its inverse, whichever is less than 



one. This topic was further investigated in some detail in Refs. |18l . Il9|, 120!, 121|, [2^. In 
particular, it was shown [3] that, by repeated and properly timed oscillator frequency 
shifts, one could reach an arbitrary degree of squeezing, even for A close to unity (neglecting 
damping, and assuming the frequency shifts to be instantaneous). A general analysis of this 
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situation, under less restrictive assumptions, was given in Ref. 

In this paper we consider the possibility of using repeated frequency shifts to produce 
GHz squeezed states in a superconducting resonant tank circuit, and the use of these states 
to improve the sensitivity of such a circuit when employed as a parametric transducer. The 



parameters of the circuit can be tuned by controlling the state of an RF SQUID inductively 
coupled to the superconducting oscillator (resonant tank circuit) 0, Is]. This distinguishes 
our proposal from Ref. [3], where the generation of squeezed states in an RF SQUID with 
a tunable junction was proposed through a single-step change of the parameters of the 
junction. 

The assumption of an instantaneous switching of the oscillator frequency is convenient 
for a proof-of-principle analysis, but not sufficient for the discussion of an experimental 
realization of the effect. Therefore, we start from the density matrix p of the tank circuit 
coupled to the superconducting oscillator, which satisfies the equation 

tdtp=[H{t),p], (3) 

where H{t) is the Hamiltonian. Additional Lindblad terms in the r.h.s. can be added to 
account for dephasing and relaxation 



Initially, H{0) = (ti;o/2)(aoao + '^o'^o)- ^ change of the oscillator frequency, Uq 



transforms the creation/annihilation operators to 171] 



a = flo H %'i ^ ~ '^0 "I '^0- 

These formulas follow from the requirement that the momentum and position operators 
do not change. This Bogoliubov transformation can be rewritten [l^] as ao ^ a = UaoU^, 
where 

U = exp[(l/4) ln(^/a;o)(4' - «o)]- (4) 
The transformation of the Hamiltonian under the unitary transformation (jlj) is given by 

H = UHoU^ - lUdtU^- (5) 

here Hq = H{0) is the initial Hamiltonian; the last term is the most important in the 



case of fast frequency changes (see, e.g., [24|). The resulting time- dependent Hamiltonian, 



expressed in terms of the original creation/annihilation operators, is 

Hit) = Ho + ^^^^^^^(«o«o + «o4 + al' + al) - ^^(4' - «o)- (6) 

Finally, by moving to the interaction representation with respect to Hq, we find the 
Hamiltonian (which, for the sake of briefness, is also denoted by H{t)) 

Hit) = "^'X^;)"" (4ao+aoat) + ^(^£ - "° (4 ^e^-o^+,2^-2..o^)_,^(42^2..ot_,2^-2..o^). 

(7) 

Hereafter, it is convenient to use the coherent state representation [3, S]. A coherent 
state I a) is the eigenvector of the annihilation operator with the (complex) eigenvalue a: 
\a): ao|a) = a\a), while (alcto = {a\a*. Each coherent state is a superposition of an infinite 
number of Fock states (states with a definite number of photons) and in the classical limit 
becomes a classical state with definite, time- dependent position and momentum (or other 
appropriate pair of canonically-conjugate classical variables). 

A density matrix p can be represented by the Wigner function W{a,a*). Examples are 
a coherent state: = |C)(CI ^ W(^{a,a*) = fe"^'""'^!^; a thermal state: Wt(«,«*) = 
^ tanh (^^^ exp — 2|ap tanh ; and a number state: p„ = \n){n\ ^ Wn{<y,a*) = 

H(zi)_g-2|apL^^4|Q,|2'j_ jjgj^g ljn{z) is the Laguerre polynomial. The complex variables a and 

4 



a* can be expressed through their real quadrature components, x^y. a = x + iy, a* = x — iy, 
which can be related to directly observable properties of an oscillator (e.g., current and 
voltage). For our purposes, the usefulness of this representation is due to the fact that there 
exists a one-to-one correspondence between the action of creation/annihilation operators 
on p and differential operations on W{a,a*): 

aop ^[a + (l/2)(9„.] W{a, a*); pao ^ [a - [1/2)8^,] W{a, a*); 
alp ^ [a* - (1/2)9J W{a, a*); pal ^ [a* + (l/2)9„] W{a, a*). 

Therefore the Liouville equation for the density matrix operator is now replaced by a partial 
differential equation for a c-number Wigner function 

idtW{a, a*) = 2^{t) [a*d^* - ad^] W{a, a*) + 2 [7(t)*a5„. - i{t)a*d^] W{a, a*) (8) 



Here we have introduced 
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l{t) 



Pit) + % 



Uj{t) 



2iuJot 



(9) 
(10) 



In the presence of dissipation, Eq. ([8]) will also contain second-order terms, and would 
be only tractable numerically. Neglecting these terms leads to an immediate simplification, 
since Eq. (|H]) is a differential equation of first order and can be solved by the method of 
characteristics. The characteristics x(t),y(t) satisfy the equations: 

1 



-X = [Im7(t)] x{t) + {f3{t) - [Rej{t)]} y{t); 



Uj = - {m + [Re7(t)]} x{t) - [Im7(t)] y{t). 



(11) 



After finding the solutions x{xo,yo,t)]y{xo,yo,t) (where xo,yo are the initial conditions) 
and inverting them to obtain xo{x,y,t);yo{x,y,t), one obtains the Wigner function for an 
arbitrary point and time from its value at t = via 



W{x + iy,x- iy, t) = W{xo{x, y, t) + iyo{x, y, t),xo{x, y, t) - iyo{x, y, t),0). 



[12) 



Eq. (fT2l) provides the complete formal solution for the quantum mechanical problem of a 
harmonic oscillator with variable frequency in terms of the characteristics of Eq. ([8]), as 



should be expected [6|, l25|]. Equations (fT2|) are of the classical Hamilton type, and can be 
solved in general only numerically. Still, a good analytical approximation can be found 
under two assumptions. First, the relative change of the oscillator frequency must be small 
(|1 — A| ^ 1). (This assumption holds for any reasonable experimental realization.) Second, 
the frequency must change either very fast (A ^ uq), or very slowly, (A ^ ujq), compared 
to the oscillator period. It was shown I22I in a general case, that only when the frequency 
is changed fast in one direction and slowly in the other, an arbitrarily strong squeezing by 
small repeated frequency shifts can take place. Therefore our analytical approximation is 
good for describing the very regime we are interested in. 

Let us first treat the fast limit. In this case, j3{t) can be neglected compared to 7(t), and 
the equations ffT^ are reduced to d = —2ij{t)a*{t); a* = 2i'y(t)*a{t). We must keep all 
terms in 7(t) until we transform this system into a second-order equation for a{t). Now, 
dropping the small terms (assuming A ^ A = f), we obtain 

a + vce — Av^a = 0. (13) 

In the case of a linear frequency change, v = const, this equation is easily solved with 
the initial conditions a{0) = a, d(0) = 2va*, yielding (see Fig. [1]) 



xoit, X, y) = Xo(t, X) = (2_,_/,)e^!.f(!^/,_2)e--. ' (14) 
yo(t, X, y) = yoit, y) = _ ^ysl^ _ ,_,,ys.v (15) 



s+ = (v^ - 1) /2; s_ = -v(Vri +1) /2. (16) 



In the slow regime, the terms with A can be neglected. The remaining terms are of the 
same order, but 7-terms are very fast (oscillating with 2ujq) and average to zero over the 
period of the oscillator. Therefore we are left with 

h = my{t): \y = -m<t). (17) 



or 



l\>^\s) - l)ds 
2 Jo 



(18) 



a = —2i(3{t)a{t), =^ a{t) = a{0) exp 

Even without solving these equations, it is clear that the slow regime cannot affect squeezing 
in any way: Eq. (|T71) describes circles in phase plane; the evolution of the Wigner function 



given by (fT2l) therefore reduces to its (slow compared to uq) rotation as a whole, without 
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changing shape. This conclusion is consistent with previous work 

The experimental realization of this proposal is rather challenging. There have been 
several reports of an ultrafast perturbation of optical microcavity modes using the dispersion 
of injected free carriers 26, 27, 3]. However, the small frequency shifts (of the order 
|A — 1| ~ 5 X 10~^) were on a picosecond timescale, which is slow with respect to the period 
of one optical cycle of the cavity modes. These processes are therefore in the adiabatic limit, 
rather than in the sudden-frequency-shift regime 17|. Thus, while they may prove useful for 
on-chip frequency conversion 29|, |30l], they are not useful for generating non-classical optical 
states. Moreover, there is little prospect for a repeated application of the perturbation 
within the sub-nanosecond lifetime of a microcavity mode. 

In optical lattices occupied by cold atoms, the relatively low oscillation frequencies (~ 1 
MHz) and precise dynamic control of the atomic potentials have allowed the demonstration 
of squeezed positional states [31[. The situation is also promising in the (O.l-lO)GHz range, 
where one can use Josephson junctions as nondissipative nonlinear elements, allowing control 
of the circuit parameters. 

If the frequency is low, so that huj — kuT, where T is 10-50 mK (dilution refrigerator 
temperatures) the so-called RF SQUID [32] configuration can be used. The system consists 
of a high-quality superconducting tank circuit, inductively coupled to a loop containing the 
Josephson contact (the single junction interferometer [s^j. In the dispersive mode of an RF 
SQUID, the effective inductance of the system becomes |32[ |: 

M2 



L 



T 



(19) 



where L^, L are respectively the self-inductance of the tank circuit and of the RF SQUID 
loop, M is their mutual inductance, and 

L 



(20) 



(3 COSV9 

is the Josephson inductance of the SQUID junction, which depends on the phase bias (p 
across it. Here (3 = 2'kLIc/^o and Iq is the critical current of the SQUID Josephson 
junction. Therefore, by varying Lp one controls the effective inductance and eigenfrequency 
of the tank: 



1 + 



^eff 



[1 + [3 cos if) 



(21) 



Here k'^ = M'^/LLt <^ 1 (in practice k ~ 0.3 can be easily realized) is the coupling coefficient 
between the tank and the SQUID loop. Since for a SQUID dispersive mode (3 < 1, the 
variations 6ujq of the eigenfrequency of the tank satisfy 6uJo ~ (0.1 — 0.01)u;o) which should 
be enough for our purposes. 

For higher frequencies hu ^ ksT, a tuneable superconducting cavity could be used. In 
such a device, a DC SQUID is incorporated to the strip resonator [33] and 6u!o ~ (0.1 — 
0.01)u;o as well. 

For both cases, the phase- changing pulse must be sharp on the scale of ujq and this is 
within the current experimental capabilities. The coherence time of the system is currently 
in the range of 1-10 /is, allowing for at least several cycles of frequency change, with the 
corresponding increase in the squeezing of the final state. 

In conclusion, we have shown that the above procedure can produce squeezed states in an 
nonlinear superconducting oscillator, so that the fluctuations of the amplitude (phase) of the 
oscillator are suppressed along certain directions in phase space, which rotate with the base 
oscillator frequency wq. By making use of this noise suppression, the measurements of the 
amplitude (phase) of these oscillators can reach a sensitivity below the standard quantum 
limit, or at least, below the thermal level. 
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FIG. 1: Squeezing of a thermal state, with ksT = 4:UJq, by repeated frequency shifts, with Amax— 1 = 
0.05. (a) Contour plots of the approximate expression of the Wigner function in the interaction 
representation (see Eqs. (12,15,17)) are shown as a function of the quadrature variables, x and 
y. Clockwise starting from the upper left: initial thermal state, the state after one, two, and ten 
cycles. The dark background corresponds to 1^ = 0, and white to W = 0.08. (b) A schematic 
frequency-versus-time dependence necessary to produce the results of Fig. la. The idle periods At 
are chosen to ensure that every fast shift occurs in the same phase with respect to the quadrature 
coordinates. Here one must account for the fact that in the Schroedinger representation the whole 
phase plane rotates around the origin with the base oscillator frequency ujq. An additional rotation 
of the Wigner function as a whole during the slow shift, Eq. (jlSp . can be neglected in comparison, 
since the additional phase \66\ < 0.5 loq t jA^^^^ — l| ^ t. 
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